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Quantum walks have been shown to have impressive transport properties compared to classical 
random walks. However, imperfections in the quantum walk algorithm can destroy any quantum 
mechanical speed-up due to Anderson localization. We numerically study the effect of static disorder 
on a quantum walk on the glued trees graph. For small disorder, we find that the dominant effect 
is a type of quantum decay, and not quantum localization. For intermediate disorder, there is 
a crossover to diffusive transport, while a localization transition is observed at large disorder, in 
agreement with Anderson localization on the Cayley tree. 
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I. INTRODUCTION 

Quantum walks constitute a promising route to the 
development of algorithms for quantum information pro- 
cessing. Successful applications of quantum walks in- 
clude a variety of search algorithms PJE], graph hitting 
problems [SHT], Boolean function evaluation [3(9], among 
others. Quantum walks are in fact a universal paradigm 
for quantum computation [10] [IT] . 

The simplest and most dramatic improvement demon- 
strated by quantum walks, when compared with the cor- 
responding classical random walks, is the hitting proba- 
bility of a walk on certain graphs. Two graphs in par- 
ticular have demonstrated an exponential separation be- 
tween the quantum and classical walks: the glued binary 
trees [5] and the hypercube [7 . An example of the former 
is shown in Fig. [I] this graph is formed by connecting 
the leaves of two binary trees of depth d. Due to the 
graph's symmetry, the quantum walk can be restricted 
to a subspace of the total (exponentially large) Hilbert 
space, known as the "column space" (also shown in Fig. 
1). It has been argued that this symmetry is the heart of 
the quantum speed-up [12]. A quantum walk algorithm 
on a slight modification of this graph, to be described 
below, was proven to have an exponential speedup over 
any classical algorithm [6]. 

Keating et al have argued that physical implemen- 
tations of these walks will not be able to maintain this 
quantum speed-up for large graphs [13]. Indeed, one ex- 
pects that physical systems will be subject to decoher- 
ence and disorder. While one can always argue that fu- 
ture quantum computers will be protected from these 
effects by error correction, it is more likely that near 
term demonstrations of quantum walks will use a phys- 
ical network. In fact, many of the recent experiments 
on quantum walks involve an encoding of the degrees of 
freedom that is not, strictly speaking, computationally 
useful [14]. Nevertheless, these experiments demonstrate 
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the dynamical speed-up characteristic of quantum walks. 
Determining when disorder or decoherence will limit the 
observability of this speed-up is an important theoretical 
problem [T5] . 

This problem may already have been encountered by 
nature. Recent evidence has shown that photosynthetic 
complexes operate using a type of quantum walk [16], 
where the interplay of decoherence with disorder plays 
a crucial role in their energy harvesting efficiency. This 
phenomenon, known as dephasing or environmentally as- 
sisted transport, has relations to earlier studies of phonon 
effects on electron transport in disordered solids [17 . Un- 
derstanding and potentially reverse-engineering this effi- 
ciency is a topic of great scientific and practical impor- 
tance. As the complexes and potential technologies under 
study can be far from the thermodynamic limit, exam- 
ining quantum transport on graphs of modest size may 
reveal new surprises. 

In this paper, we examine the continuous-time quan- 
tum walk on the glued trees graph with static disorder, 
corresponding to the usual Anderson model with dis- 
ordered on-site energies [IB] . Previous authors [13] in- 
troduced disorder through an effective one-dimensional 
representation (the column space). Using the fact that 
all eigenstates are localized in one-dimension with arbi- 
trary disorder, they concluded that Anderson localiza- 
tion would cause an exponential suppression of the hit- 
ting probability. From a physical perspective, however, 
this model has some flaws. First, their model introduces 
a highly correlated form of disorder, in that all of the 
on-site energies on a given column are the same Second, 
the glued trees graph is quite similar to a Cayley tree, 
whose dimensionality is formally infinite, as far as An- 
derson localization is concerned [19]. The Cayley tree 
exhibits a localization transition at large disorder, and 
thus significant speed-up may still be possible for weak 
disorder. 

We have performed a numerical study of this prob- 
lem for the full Hilbert space of the glued trees graph 
and modifications thereof. Analysis of the eigenvalues 
and eigenvectors agree with previous studies of localiza- 
tion the Cayley tree. We pay special attention to the 
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less well studied regime of weak disorder, using dynam- 
ical simulations to find a type of quantum decay of the 
quantum walk from the column space. The resulting hit- 
ting probability is well simulated by a model with local 
decay from each column state. We further consider the 
crossover from wavelike to diffusive transport for inter- 
mediate levels of disorder. A scattering theory analysis 
of this regime indicates a transition from the quantum 
walk to a classical random walk. 

Our results augment the many results for one- 
dimensional quantum walks with disorder. These include 
theoretical studies for the continuous-time quantum walk 
[20] and the discrete-time quantum walk [21], and the 
many recent experiments on atoms [22], ions [23j [24] . 
and photons [25j [26] . Higher-dimensional quantum walks 
could be realized in these or other systems, such as net- 
works of superconducting circuits [27] [28]. In particu- 
lar, our simulations show that the effects of disorder on 
quantum transport can be identified in systems of 10-100 
lattice sites. 

This paper is organized as follows. In Section II we 
review the known results for Anderson localization on 
the Cay ley tree, and in Section III numerically study the 
phase diagram for the two types of glued trees, and find 
evidence of the localization transition. In Section IV we 
study the time-dependence of the quantum walk with dis- 
order, and introduce the local decay model for weak dis- 
order. In Section V we consider a scattering theory model 
for transport through the glued trees graph, and find ev- 
idence for the transition to classical random walk. We 
conclude in Section VI, comparing our results with the 
hypercube and highlighting the major open questions. 
Results for the quantum scattering and (classical) diffu- 
sive transport through the glued trees are found in the 
Appendices. 



II. LOCALIZATION AND DIFFUSION ON THE 
BINARY TREE 

Anderson localization is the phenomenon that, for a 
sufficiently large amount of disorder, the eigenstates of 
a quantum system become exponentially localized about 
the nodes in a graph [18 . Localization transitions are 
a rich field of study [19], for which both symmetry and 
dimensionality play key roles. The quantum walk we 
consider corresponds to the simplest such model, a tight- 
binding Hamiltonian with random on-site energies 

where 7 is a hopping rate, ej is a set of random variables, 
uniformly distributed in the range [— W/2, W/2], and ct 
is the creation operator for an excitation at site j, and 
the sum is over neighboring sites. This model, originally 
inspired by electron transport in solids, can be used for 
many physical systems, such as quantum state transfer 
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FIG. 1: Glued binary trees graph G4, showing the reduction 
to the column space for both the classical (top) and quantum 
(bottom) random walks. 



of a single excitation on a qubit network [27 , 29, 30 or 
exciton transfer in a photosynthetic complex [T6] . 

Much is known about the Anderson model given by 
Eq. [l] [19 . For example, for a one-dimensional infinite 
lattice, localization occurs for arbitrarily small amounts 
of disorder [31] . It was argued that this property was 
generic to quantum walks [13]. However, it is known that 
for systems with dimension greater than two, localization 
requires a sufficiently large amount of disorder, i.e. there 
is a localization transition as the disorder is increased 
[32] . At first glance, the glued trees graph of Fig. 1 
might appear to be a subset of a two-dimensional system, 
and thus one might expect localization for even small 
amounts of disorder. 

In fact, the infinite binary tree, or more generally the 
infinite Cayley tree (also known as the Bethe lattice [33 ) 
has been used as a model for an infinite-dimensional sys- 
tem. This has been studied extensively, and exhibits a 
localization transition found numerically [17, 34-36, and 
by an analytical mean field calculation [37j [38] . For this 
system, there is a mobility edge in the energy spectrum, 
such that eigenstates inside the mobility edge are ex- 
tended, while those outside are localized. For sufficiently 
large values of disorder, there is a localization transition 
beyond which all eigenstates are localized. For a binary 
tree (or Cayley tree with K = 1), this transition occurs 
for W c « 17 QI1 MM EH] . 

The description of localization given above, in which 
the eigenstates of the system exhibit a transition from 
extended to localized states, can be called eigenstate lo- 
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calization. There are two other ways to identify local- 
ization that we will encounter in this paper. The first 
is dynamical localization, in which the spreading of a 
wavepacket shows a saturation as a function of time. The 
second is spectral localization, in which the eigenvalues of 
the system move from an absolutely continuous spectrum 
(corresponding to extended states) to a pure point spec- 
trum (corresponding to localized states). These other 
indicators of localization, which have also been studied 
extensively, we now briefly summarize. 

The spectral properties of the system were studied nu- 
merically [36j [39], and found to exhibit a transition in 
agreement with the self-consistent approaches described 
above. These numerical studies are sensitive to the han- 
dling of the boundary of the tree [40], a fact to which 
we will return in Section III. While there are some sub- 
tleties regarding the the phase diagram at weak disorder 
[4lJ[42], the existence of a localization transition for large 
disorder is well established. For small disorder, the ex- 
istence of extended eigenstates has been proven [43]. It 
has also been proven that, for small disorder, states that 
are initially localized spread ballistically [44] . 

The ballistic spreading does not preclude classical be- 
havior, however. The random walk on the Bet he lattice 
has been studied, and can be mapped onto an asymmet- 
ric random walk on a one-dimensional half-line [45], as 
indicated in Fig. [I] (e.g. the left half). A classical walker 
is twice as likely to move right as left, and this asymme- 
try leads to the peculiar fact that the classical walk also 
spreads ballistically. Exact results and limits have been 
established for this and other properties of the classical 
walk [46[ |47] . Note that the asymmetry is linked to the 
exponential growth of sites away from the origin, such 
that one often calls the Bethe lattice a graph of infinite 
dimensionality. 

III. QUANTUM WALK EIGENSTATES 

The continuous-time quantum walk [48j [49] is precisely 
the quantum dynamics of a single particle moving on a 
graph. This is given by the Schrodinger equation 

.d4_l = _^ A . k ^ (2 ) 

k 

where Ajk is the adjacency matrix for the graph and 7 is 
hopping rate. One could also use the Laplacian matrix, 
or introduce potentials to implement search algorithms 
[2j[4], but here we consider the addition of static disorder, 
such that 

i^jr = -l^Ajk^k^e^j, (3) 

k 

where the on-site energies €j are i.i.d. random variables 
uniformly distributed in the range [— W/2, W/2]. 

In this section we consider the nature of the eigenstates 
of the Hamiltonian of the quantum walk with disorder, 



namely H\V) = E\V), with U = H + W, with the 
unperturbed Hamiltonian given by 

N d 

H = -7 E A Jk\j)(kl (4) 
and diagonal static disorder 

N d 
3 = 1 

where Nd — 3 x 2 d — 2 is the number of vertices for the 
glued trees graph G^. We begin by analyzing Ho, and 
follow by studying the eigenvalues and eigenstates of H. 

A. Eigenstates without disorder 

As described above, there is a great deal of symmetry 
in the glued trees graph G^, and a great deal of structure 
in the eigenstates and eigenvalues of the system. We be- 
gin by providing a notation for the graph. We consider a 
labeling along the "columns" and "rows" of the graph, of 
the form (j, n), where j = — >• 2d indicates the distance 
from the left root along the graph, and n = —> Nj^ — 1 
is the location within a given column. Here Nj^ is the 
number of sites in a given column j, given by Nj^ — 2 J 
for j < d, and Nj^ = 2 2d ~ J for j > d. The coordinates 
( j, n) can be combined into a single coordinate v by the 
following rule 

_ / 2? + n for < j < d, , , 

V \ 3 x 2 d - 2 2d+1 -i + n for d < j < 2d. W 

Note that v ranges from 1 — » 3 x 2 d — 2. 

The adjacency matrix elements A VyW are equal to one 
if vertices v and w are connected, and zero otherwise. 
This can be given in terms of the coordinates ( j, n) by 
observing that (j, n) is connected to 

(j - 1, [n/2j ), (j + 1, 2n), (j + 1, 2n + 1) for j < d, (7) 

and 

(j + 1, [n/2j ), (j - 1, 2n), (j - 1, 2n + 1) for j > d. (8) 

What is most important is that a given vertex is sym- 
metrically coupled to those sites one step further from 
the left root (or closer to the right root). This allows 
us to express the eigenstates of the system in terms of 
"column-states" that are equal superpositions of states 
on a given column j. This column-space reduction is 
well-known for its utility in the analysis of the quantum 
walk [5] El [[2]. 

Letting |j, n) denote the Hilbert-space vector associ- 
ated with vertex (j, n), we define the column-space vector 
I col j) by 

|colj> = -== l^>- W 

V iV ^ n=0 
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By the symmetry noted above, the Hamiltonian Hq = 
A acts on the column-space states as 

Ho|col j) = -V^Tlcol j-l)- V^Tlcol j + 1), (10) 

where one factor of y/2 is due to the number of connec- 
tions to a neighboring column, and the other due to the 
normalization of the column states. Hence, we can re- 
duce the dynamics to a quantum walk on a finite line, 
whose eigenstates are equally well-known 



1 



2d 
3=0 



sin 



k (j + 1) 7T 

2(d+l) 



with energies 

Ekd = — 2V27 cos 



kit 



2(d+l) 



Icolj), (11) 



(12) 



and k = 1 2d + 1. We have annotated the states by 
the depth of the graph, as there are in fact many more 
eigenstates for G^, whose enumeration we now consider. 

The glued trees graph is self-similar, in that Gd con- 
tains 2 V subgraphs, each equivalent to Gd-v These sub- 
graphs can be grouped into 2 V ~ X pairs, each pair formed 
by removing the roots of a larger graph equivalent to 
Gd-v+i- That is, we can repeatedly split the glued trees 
graph by removing the left and right roots, such that 
Gd contains 2 copies of Gd-i, 4 copies of G^-2, and so 
forth, formed by removing the roots until we are left with 
2 d copies of Go (the isolated vertices at the center of 
the graph). Some representative subgraphs of G4 are 
shown in Fig. [2] On their own, each subgraph would 
have eigenstates of the form of Eq. (11). To find how 
these subgraphs contribute to the spectrum of G^, we 
observe that an equal but opposite-signed superposition 
over two paired subgraph eigenstates is an eigenstate of 
Gd- This occurs because the components with opposite 
phase, when acted upon by Ho, will interfere destruc- 
tively on the two nodes to which they are connected on 
the larger graph. By this pairing, we can thus construct 
the complete set of eigenstates for G^. 

To see this more clearly, we define a set of "sub- 
column" states whose elements combine the paired sub- 
graphs described above. These are given by 



|scolj;a,i/) 



(2a+l)AT d _ 

£ 

n=2aN d - 



\j + v, n) - \j + v,n + N j4 - U ) 

V 2N 3,d-v 



(13) 

where j = — » 2(d — v) indicates the column in Gd-v> 
a = — » 2 v ~ l — 1 labels the paired subgraphs, and v = 
1 — >• d indicates the depth of the subgraphs' left root. 
Using these states, the remaining eigenstates of Gd can 
be obtained by using the eigenstates \^k,d-u) from Eq. 
( fTTj ) with I col j) replaced by |scol j;a,v) and eigenvalues 
E k ^ d -v from Eq. ( fl2| ) ? where k = 1 -» 2(d - v) + 1. 
Defining 




FIG. 2: (Color online) Glued binary trees graph G4, with 
several highlighted subgraphs G3 (top, in gray) and two copies 
of G2 (middle, in red and bottom, in blue). 



the total spectrum can then be written as 



2 „- 



(15) 



crd = {E k , d ,k = 1 -> 2d}, 



(14) 



This spectrum exhibits a very large degeneracy, espe- 
cially for E = 0, which has a multiplicity of 2 d (one from 
each copy of <Jd-v> and half from the cr ). This large de- 
generacy (also observed in [40]) occurs for trees that are 
both glued and unglued and can make numerical analysis 
of the disordered system problematic, as will be described 
below. 



B. Eigenstates with disorder 

The introduction of static disorder changes both the 
eigenvalues and eigenvectors of the system. For the Cay- 
ley tree, early studies used a self-consistent approach 
[3U|38] to find the localization transition and a phase dia- 
gram between extended and localized eigenstates. Recent 
numerical studies [35] [36] have confirmed these earlier re- 
sults, which we now summarize. 

The eigenvalue spectrum has been studied for a single 
Cayley tree numerically by [39] through the use of spec- 
tral statistics. A transition of the distribution of the en- 
ergy level spacings from Wigner to Poisson (indicative of 
a localization transition) was observed when the tree was 
modified so that the leaves are randomly connected to 
each other. This random connection presumably reduces 
the prevalence of the zero eigenvalue described above, 
which would otherwise lead to a Poisson distribution for 
the spectral statistics [40] . We performed a similar anal- 
ysis for the glued trees graph, with the leaves connected 
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as in Fig. [T] which we call a simple glued trees (SGT) 
graph, or interrupted by a random cycle (as in [6]), which 
we call a modified glued trees (MGT) graph. An exam- 
ple of the MGT (with a regular connection [50 ) is shown 
in Fig. 9. A Wigner to Poisson transition was observed 
for the MGT, while the SGT exhibited a spectrum that 
was always far from Wigner. For the remainder of this 
section, our results were obtained using the MGT. 

The localization of the eigenstates and the localization 
transition can be found by studying a simple property 
of the eigenstates, namely the inverse participation ratio 
(IPR) h- 



(16) 



where we assume that the states are normalized with 
h(^) — l^jl 2 = 1- This quantity has the property 
that an eigenstate localized to a single site has h(^) = 1, 
while an eigenstate extended over N sites has I2WO = 
1/N. We further define an averaged value of this quantity 
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N(E;AE) 



W), (17) 



\(^\nW-E\<AE 



where N(E; AE) is the number of eigenstates found to 
have eigenvalue Ej with in the range E — AE < Ej < 
E + AE, leaving the dependence on AE implicit.. By cal- 
culating hiE) the eigenstates of the system as a function 
of energy and disorder, the phase diagram and localiza- 
tion transition can be visualized. 

For the SGT, we again encounter a difficulty in that 
there are a large number of states with an h^) = 1/2, 
associated with the 2 d ~ 1 states in cr (with E = 0). How- 
ever, by using the MGT, we can calculate meaningful 
results for 12(E) for various disorder strengths W to ob- 
tain the diagram in [3] Here we have let d = 8 and set 
AE = 0.157 and averaged over 500 realizations of H. 
We observe a gradual movement of extended states (with 
small I2) to localized states (with I2 > 1/2) as disorder 
is increased. Also shown are the expected results for an 
infinite Bethe lattice, obtained using the self-consistent 
method of Miller and Derrida [34] 

To obtain an estimate of the localization transition, 
we fix our attention to states near E = 0, and repeat 
the calculation of the averaged IPR for graphs of various 
sizes. As expected, ^(0) exhibits a small size-dependent 
value for W = 0, which increases to approximately 0.5 for 
large W. At a certain value, the averaged IPRs for all of 
the graphs coalesce, from which we estimate W c « 17, in 
agreement with recent results for the Cayley tree [35| f36] . 



IV. QUANTUM WALK DYNAMICS 

The dynamics of a quantum walk on the glued trees 
graph with disorder has been studied using the one- 
dimensional column-space model by Keating et al. [13]. 
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FIG. 3: Inverse participation ratio 12(E) as a function of 
energy and disorder, for d = 8. For each value of disorder 
the IPR was averaged over 500 realizations of 7i with AE — 
0.15 (with 7=1). Also shown is the mobility edge from the 
self-consistent theory, predicting a localization transition with 
W c « 17. 
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FIG. 4: (Color online) Inverse participation ratio 12(E) at 
the band center (E — 0) as a function of disorder for various 
depths. For each value of disorder the IPR was averaged over 
the middle 100 eigenvalues and a number of realizations given 
by 500, 250, 125, and 50 for d = 5, 6, 7, and 8, respectively. The 
localization transition occurs when the curves collapse near 
W « 17. 



In the previous section, however, we have seen that the 
eigenstates of the MGT undergo a localization transition 
at large disorder. This leaves open the possibility that a 
dynamical speed-up can be observed for small disorder. 
To explore this possibility, we again turn to numerical 
studies. There are three issues to be studied: first, what 
is the probability to hit the right root should the quan- 
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turn walk begin at the left root? Second, how does this 
probability decay as a function of disorder and size, and 
by what mechanism? Finally, if the disorder quantum 
walk does not hit the right root, how far does it get? 

The first question can be answered by calculating the 
average of 



Pbit(t) = \(col 2d\e~ mt \co\ 0)| 



(18) 



for instances of the SGT Hamiltonian with disorder 
strength W and for various sizes d. This corresponds to 
the single-shot measurement procedure [7 for the quan- 
tum walk, and a representative set of Phit(^) as a function 
of time and values of W are shown in Fig. [5j The first 
thing to observe is the oscillatory structure, due to the 
Bessel function structure for £>hit (t) IMS , from which 
we can see that probability is maximized at the hitting 
time 



thit 



2d + 1 + 1.0188 d 



IX 1/3 

V 



(19) 



with a value pd ~ d~ 2 / 3 for large d. Second, as disor- 
der increases the maximum of the hitting probability is 
seen to decrease. This hitting probability can be approx- 
imated by 




FIG. 5: (Color online) Hitting probability PhitM as a func- 
tion of time for the quantum walk on the glued trees graph 
with d = 15 and various values of disorder. The symbols were 
calculated using exact numerical simulations with (from top 
to bottom) W = 0,0.2,0.4,0.6,0.8, 1.0, and 1.2, with 10 re- 
alizations of % for each value of W . The solid curves were 
calculated using the locay decay model (see text). The inset 
shows the hitting probability phit(£hit)/P<i as a function of d 
with symbols from numeri cal simulations and the curve from 
the approximation of Eq. (20) (see text). 



Phit(thit) -Pdexp 



(20) 



The exponential decay of Phit/Pd is shown in the inset to 
Fig. 

By what mechanism does this decay occur? It is not 
associated with eigenstate localization, as the eigenstates 
of the system are well within the extended regime in the 
phase diagram of Fig. [3J To understand this, we turn to 
another quantity of interest, the column-space probabil- 
ity 



2d 



(21) 



Quantum walks on the glued trees graph that are initially 
in a column space state |col jo) will remain in the column 
space in the absence of disorder. Once static disorder is 
introduced, the eigenstates lying within the column space 
become coupled to the other eigenstates of Gd- These 
are associated with the subgraphs of G^, and have zero 
amplitude on the graph's two roots. The resulting decay 
of p co \(t), shown in Fig. pleads to the decay of phit(t)- 

We can provide an analytical estimate of this decay 
through perturbation theory. Assuming the walk begins 
in a column state | col jo ), taking the second order ex- 



pansion of p co \(t) gives 

2d 

Pco\(t) = ^2 I (col j |exp(-z H t) I col j ) I 2 

3=0 

2d 

« £ |{ col j I col jo >| 2 + * 2 |{ col i I ^ I col io )| 2 

3=0 

-t 2 « col j I col jo ) ( col j I H 2 I col j )) , 

(22) 

where the Hamiltonian is given by H = Hq+H'. It is 
straightforward to show that Ho contributes nothing to 
the quadratic term, and for %' given by Eq. ^ we find 



PcolCO = 1 



N 2 , 

Jo,d 



i i^j 



Averaging p co \ (t) over the disorder we find 



(23) 
(24) 



where we have used (e^) = ^W 2 (5^. 

This result for the short time decay applies to any 
graph for which the quantum walk can be projected onto 
a column space. On longer timescales, we conjecture that 
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the decay will have an exponential character, while re- 
taining the position dependence. Hence, extrapolating 
from the short time result, we postulate a model of "lo- 
cal (exponential) decay", in which the walk evolution is 
computed as in the ideal column space representation, 
but the probability at each site j is allowed to decay: 



Pj(t) =Poexp 



t W 2 
12~ 



1 



1 



N, 



3,d 



This is equivalent to applying the mapping 



H H 



col 



■iT/2, 



(25) 



(26) 



where H co i is the projection of the graph Hamiltonian 
onto the column space, and T is given by 



r = 



12~ 



2d 

£ 

3=0 



|colj)<colj|, (27) 



with 7 denoting the unit time hopping probability from 
Hq. This expression for T can be anticipated by applying 
Fermi's Golden Rule to W in the eigenstate basis, but the 
presence of a large discrete component to the spectrum 
(<Jo) prevents us from making a precise derivation. Nu- 
merical evidence, however, suggests that this is in fact 
the correct mechanism. 

Along with the exact numerical results for phit an d p C o\ 
in Figs. |5]and|6j we have included results from the local- 
decay model calculated in the (exponentially smaller) 




FIG. 6: (Color online) Column space probability p co \(t) as 
a function of time for the quantum walk on the glued trees 
graph with d = 15 and various values of disorder. The sym- 
bols were calculated using exact numerical simulations with 
(from top to bottom) W = 0, 0.2, 0.4, 0.6, 0.8, 1.0, and 1.2, 
with 10 realizations of H for each value of W . The error 
bars indicate the standard deviation of each average. The 
solid curves were calculated using the local decay model (see 
text). 



column-space representation for Hq. We observe that 
the local decay model predictions (solid lines) agree well 
with the results of simulations (points), to within one 
standard deviation. Our model accurately reproduces 
many features of the exact simulations, such as the vari- 
ation in the decay rate of p co \ and the oscillations in phit • 
Furthermore, if we use the column space probability to 
predict the hitting probability at the target node (oppo- 
site root), the agreement is quite good. This supports 
the claim that the disorder-induced reduction in quan- 
tum transport is primarily explained by decay from the 
column space. 

So far we have looked at the probability at the right 
root and the column space, but a more global characteri- 
zation of the walk propagation can be found by analyzing 
the average depth reached by the quantum walk, 



r(t) = (m\r\m), 
where the column position operator f is defined as 



2d N jtd -1 
j=0 n=0 



n\ 



(28) 



(29) 



with the property r\ col j) = j\co\ j). The value of 
r(t) thus gives a snapshot of the expected position of the 
walk along the graph at any one time. This is displayed 
in Fig. [7] as a function of time and disorder. In the 
absence of disorder, the average depth shows an oscilla- 
tory character consistent with ballistic propagation of the 
wavepacket and reflections at the two ends of the graph. 
Disorder-induced decay from the column space causes the 
amplitude of the oscillations to decay faster than in the 
ideal case, ultimately causing a "damping" of the oscil- 
lations. Therefore, the walk has a reduced probability 
of traversing the graph, and substantial probability is 
instead deposited in the center of the graph, where the 
concentration of nodes is highest. 

In Fig. IV we plot the maximum value of r(t) in the 



range t < 3£hit against the strength of disorder, illustrat- 
ing the localization transition on the glued trees graph. 
We see that for small disorder (W < 2), the maximum 
value of r is high, corresponding to the quantum walk 
hitting the right root, for which r = 2d. As disorder 
increases, the curve falls, and begins to level off at the 
graph center (r « d for 2 < W < 4). For larger amounts 
of disorder (W > 4), the curves continue to decrease, 
converging on a single value around W = 16, precisely 
when we expect all eigenstates to be localized. 

This analysis of the quantum walk dynamics strongly 
suggests a type of quantum-to-classical or "wavelike-to- 
diffusive" crossover at weak disorder. Note that this not 
a "ballistic-to-diffusive" crossover, and does not conflict 
the ballistic spreading of wavepackets found by Klein [44] 
for the Bet he lattice at weak disorder, as classical diffu- 
sion also leads to ballistic spreading (r(t) ~ t). It does, 
however, suggest that there may be a length scale (the 
mean-free-path) that limits the size of graphs for which 
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FIG. 7: (Color online) Average distance r(t) for d = 15 as a 
function of time and disorder W. For each W, the simulations 
were averaged over 10 realizations of 7i. The quantum oscilla- 
tions decay in time for small disorder, with a critical damping 
near W ~ 2, indicating a type of quantum-to-classical transi- 
tion. 




FIG. 8: (Color online) Maximum average distance r(t) as a 
function of disorder W for various depths. For each value of 
disorder the maximum of r(t) was averaged over 1000, 100, 
and 10 realizations for d = 5, 10, and 15, respectively. For 
intermediate disorder (2 < W < 4), there is a quantum- 
to-classical transition, while for large disorder (W > 15) a 
localization transition is observed. 



a speedup could occur in the presence of disorder. That 
is, the exponential decay of the hitting probability (in d) 
seen in Fig. [5] may be interpreted as the classical proba- 
bility for a walker to traverse d sites given a mean-free- 
path of order i ~ 1/W 2 . To understand this crossover in 
more detail we extend our analysis of the quantum walk 
to a transport model. 



Incident Wave 



Transmitted Wave 




Reflected Wave 



FIG. 9: Scattering approach to the quantum walk, in which 
an incident wavepacket is transmitted and reflected along the 
"tails" connected to the modified glued trees graph. 



V. QUANTUM WALK TRANSPORT 

To identify the quantum-to-classical transition, we 
consider the quantum transmission through the glued 
trees graph, subject to disorder. To transform the quan- 
tum walk into a transmission problem, we attach "tails" 
to the input and output nodes, and look at the transmis- 
sion coefficient through the graph for a wavefunction of 
the form 



gi/en _|_ 
'~j~^ikn 



Ke- ikn for n < 

for n > 2d - 



1 



(30) 



This represents an ingoing wave that is reflected and 
transmitted through the graph, as illustrated in Fig. [9] 
This type of quantum walk was been used to develop a 
quantum algorithm for NAND-tree evaluation [8 , and 
has been generally analyzed in [51] . We consider the 
transmission probability T — |T| 2 as a function of depth 
and disorder. 

Using a standard analysis for transmission in tight- 
binding lattices [52], we find the transmission amplitude 



T = (col 0|-, 



2i sin k 



H — 2 cos k 



| col 2d+ 1) 



(31) 



where 



U = H + e lk (|col 0)(col 0| + |col 2d + l)(col 2d + 1|) 

(32) 

is an effective Hamiltonian for the MGT graph alone 
(note that this graph has 2d + 1 columns). This quan- 
tity can be calculated by diagonalizing the non-Hermitian 
Hamiltonian and forming the appropriate matrix ele- 
ments in T. 

The resulting transmission probability T = |T| 2 is 
shown as a function of momentum k and disorder W in 
Figs. 10 and [II] for depths d = 5 and d = 6, respectively. 
For small disorder, there is a size-dependent oscillatory 
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t = in 2 




FIG. 10: Transmission probability T as a function of momen- 
tum k and disorder W for a modified glued trees graph of 
depth d = 5. For each value of momentum and disorder, the 
transmission was averaged over 250 realizations of 7-L. 




FIG. 12: (Color online) Transmission probability T for k = 
tv/2 as function of disorder and for various depths. The var- 
ious symbols are d — 7 (blue circles), d — 8 (red upward 
triangles), d — 9 (blue squares), and d — 10 (red downward 
triangles), averaged over 500, 200, 100, and 50 realizations of 
respectively. The solid black curve is the transmission 
probability for the classical random walk T c (see text). 



t = \r\ 2 




der can be understood using a classical model, described 
in the Appendix B. This model uses a diffusion constant 
proportional to the me an- free- path A ~ £ ~ W~ 2 and 
leads to a classical transmission probability of the form 



To 



i + c{whf 



(33) 



where the coefficients T and c presumably depend on 
the exact mapping of the disordered quantum walk to 
the diffusion equation, such as the method of [53 or the 
results of [54 . This expression for Eq. (33), fit using 
Tq = 0.8 and c = 0.2, is shown in Fig. 12 This agree- 
ment, for intermediate values of disorder, provides con- 
firmation of the quantum-to-classical crossover observed 
in the dynamical studies of the previous section. 



FIG. 11: Transmission probability T as a function of momen- 
tum k and disorder W for a modified glued trees graph of 
depth d = 6. For each value of momentum and disorder, the 
transmission was averaged over 100 realizations of %. 



structure as a function of k due to resonances, much like 
those found in [52] (and briefly described in Appendix 
A). These oscillations in T disappear when the disorder 
strength W ~ 2, after which the transmission decays 

shows T for k = tt/2 as a func- 
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monotonically. Figure 
tion of disorder for many graph sizes, which all exhibit 
the same behavior for W > 2, indicative of a transition 
in T. 

The decay of the transmission probability with disor- 



VI. CONCLUSION 

In this paper, we have carried out an investigation of 
the effects of diagonal disorder on quantum walks on the 
glued trees graph. While disorder does lead to localiza- 
tion in the strong disorder limit, we find the primary 
effect in the case of small disorder to be quantum decay 
out of the column space. The quantum decay can be ac- 
curately modeled in the column space by a non-unitary 
mapping that enforces position dependent decay of the 
probability. This local decay model is efficient to com- 
pute, owing to the exponential reduction in size of the 
representation, yet it allows prediction of the end-to-end 
hitting probability, and should be extendable to other 
graphs with similar symmetries. 
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One such graph is the hypercube. This problem had 
been previously studied for quantum state transfer [27] , 
where the effects of off-diagonal disorder were empha- 
sized. Numerical simulations for diagonal disorder, how- 
ever, provide very similar results to those found in Sec. 
IV, with one important difference. The hitting time for 
the hypercube is independent of the dimension d (which 
is analogous to the depth of the glued binary trees graph). 
The local decay model then predicts that the hitting 
probability should decay as e~ w Z 12 , a result borne out 
by simulations. Thus, for the hypercube, quantum trans- 
port outperforms classical transport for small W. 

Such a result is also possible for the glued trees graph. 
The exponential suppression of the hitting probability 
occurs for the specific case of a state initially localized to 
the left root of the graph. By using a graph with "tails" , 
the results of Sec. V show that appreciable transport 
is possible for W < 2, provided there is no disorder in 
the tails and an appropriate initial state can be found. 
In addition, alternative measurement strategies [55], or 
the inclusion of traps [49] could provide opportunities 
for speedup. Exploring this possibility could provide ad- 
ditional context for understanding environmentally as- 
sisted quantum transport [T6] . 

In summary, we have performed an analysis of the ef- 
fect of static disorder on a quantum walk on the glued 
trees graph. For small disorder, we find that the dom- 
inant effect is a type of quantum decay, and not quan- 
tum localization. For intermediate disorder, there is a 
crossover to diffusive transport, while a localization tran- 
sition is observed at large disorder, in agreement with 
Anderson localization on the Cayley tree. Our results 
suggest that intermediate disorder will inhibit any quan- 
tum speedup, but also that large speedups are possible 
for quantum walks on complex networks with small dis- 
order. 
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Appendix A: Quantum Walk Transport 

The transmission through the modified glued trees 
graph can be calculated by using the following ansatz 
for the column-space wavefunction 



MGT 



(») 



gifen _|_ — ikn 

Ae ikn + Be- fkn 

Ce ikn + J) e -i~kn 
'j-'gikn 



for n < 0, 
for < n < d, 
for d + 1 < n < 2d - 
for n > 2d+ 1. 



1, 



(Al) 

To determine T, one must use continuity of \I/ at n = 
and 2<i+l, and the Schrodinger equation at n = 0, d, 



and 2d + 1, with E = —2jcosk = —2^7 cos /c. These 
provide six equations for the six unknowns A, B, C, D, 7Z, 
and T. For general k, these equations are most conve- 
niently solved by computer, and have a rather compli- 
cated solution. For k = tt/2, however, the solution for T 
is remarkably simple 



T(k = tt/2) = 



9 + (-l) d " 



(A2) 



This oscillation of T(k = tt/2) as a function of d for the 
MGT facilitated the classical interpretation of Fig. 12. 

For comparison, the appropriate ansatz for the simple 
glued trees graph is 



SGT 



(») 




for n < 0, 

for < n < 2d, 

for n > 2d, 



(A3) 



and the corresponding solution for the transmission co- 
efficient is T{k = tt/2) = 1, independent of d. 



Appendix B: Classical Walk Transport 

The scattering theory approach of Sec. IV can be mod- 
ified to study classical diffusive transport through the 
modified glued trees graph. The model, depicted in Fig. 
[13} consists of an ingoing particle flux T to the left root. 
These particles then undergo a random walk in the inte- 
rior of the graph, with diffusion rate A, but can escape 
from either the left root (with rate Ai e ft) or the right root 
(with rate A r i g ht)- This is indicated by the directed paths 
on the left and right "tails" on the graph. In the steady- 
state, there is an outgoing flux to the left ri e f t = Pieft^ieft 
and an outgoing flux to the right r righ t = PrightA r i g ht, 
with T = ri e ft + Tright • The classical transmission coeffi- 
cient introduced in the text is T c = r r i g ht/r. 

The calculation of T c can be found by solving for the 
steady-state of the master equation 



dpj_ 
dt 



— 5^ Lj^Pk 

k 



(Bl) 



where Lj t k is a matrix of rate constants, and bj is zero 
except for the input site (where it equals T). Due to 
the symmetry of the modified glued trees graph, we can 
write this equation in terms of the probabilities for the 
walker to be on each column, so that j = — >• 2d + 1, 
Pieft = Po, bright = V2d, and bj = T5 jj0 . One can find the 
steady-state probability distribution from p ss = —L~ x b. 

Here we provide an analytical calculation by solving 
for the steady state by recursion relations. In full, the 
master equation reads 



dpp 
dt 



-(2A + A left )po + Api+r, 



(B2) 
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Requiring these two to solve dpd/dt = 0, we obtain one 
equation for po and P2d+i- 



2 d + (2 d -l)^f, 



Po 



2 d + (2 d - 2~ 1 )^± 



(2 d -l)f 

This can be combined with the conservation law 

AleftPO + A rig htP2d+l = T, 



(B9) 
(BIO) 



FIG. 13: Transport approach to classical random walk on the 
modified glued trees graph. A steady-state flux is F intro- 
duced to the left root and carried away by the "tails" on the 
left and the right, with F — F\ e f t + r rig ht- 



d i 2Ap J _i - 3Xpj + Ap i+ i for < j < d 
3 — ' 2\pj_i — 3Xpj + 2Ap J+ i for j = d and d + 1 



dt 



and 



Xpj-i — 3Xpj + 2Ap J+ i for d + 1 < j < 2d + 1 

(B3) 



dp2d+i 
dt 



Xp 2 d ~ (2 A + Aright )jP2d+i- (B4) 



For the steady-state dp/dt = 0, we can solve the rec- 
curence relations 

Pj+i = 3pj — 2pj-i for j = 1 — >> d (B5) 

Pj-i = 3pj — 2p J+ i for j = 2d + 1 — » d (B6) 

in terms of po an d P2d+i? respectively. We find 

Pj = 2*po + - l)^po - - 1)^ (B7) 

to-i = 2> 2i+1 +(?-l)^ te+1 (B8) 



to solve for p and P2d+i- For the latter, we find 



P2d+1 



r/A 



left 



1+ (^S i ) +2 ( 1 - 3x2 " <i " 2 )(^ i ) 



Forming T c = A r i g htP2d+i/r, and simplifying we find 



(Bll) 



1+ (5i) +2 ( 1 - 3x2 " d " 2 )(^ & ) 



(B12) 



To complete our calculation, we estimate the diffusion 
rate A as the product of the hopping rate (7) and the 
transition probability between sites (7 2 /W^, for large 
W). Finally, setting Ai eft = A righ t = 7 and A = 7 s / W 2 
yields 



1/2 



l + (l-3x2- d - 2 )(W/ 7 ) 2 ' 



(B13) 



a result very similar to that used in Sec. VI. 
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